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ABSTRACT 

The  influence  of  internal  heating  by  mesoscale  eddies  on  the 
oceanic  temperature  structure  is  tested  by  including  a  vertical  heat 


transport  term,  w'T',  in  the  ocean  thermal  budget  equation.   The 
value  for  this  term  was  determined  from  a  theoretical  study  by  Gill, 
Green  and  Simmons  (1974).   In  the  past  studies  of  the  thermocline 
problem,  the  heat  transport  term  has  been  neglected.   Analytical  and 
numerical  test  cases  were  run  solving  the  heat  equation  with  and 
without  the  mesoscale  eddy  heat  transport  term  for  baroclinic  and 
barotropic  conditions  in  the  ocean.   Whenever  heating  by  mesoscale 
eddies  was  included  in  the  heat  equation,  excessive  cooling  occurred 
in  the  upper  1000  m  of  the  ocean.   The  results  of  the  test  cases  in- 
dicate that  the  eddies  do  have  a  distinct  influence  on  the  ocean's 
thermal  structure  and  need  to  be  included  in  solutions  to  the  heat 
equation  for  the  ocean. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 9 

II.  EQUATION  AND  MODEL  DEVELOPMENT  --------------  13 

A.  HEAT  EQUATION 13 

B.  M  EQUATION I4 

C.  BOUNDARY  CONDITIONS- 16 

D.  EDDY  HEAT  TRANSFER -_- _ 18 

III.  ANALYTIC  AND  NUMERICAL  RESULTS  -------  --  -  23 

A.  VERTICAL  MOTION  AT  THE  SURFACE  EQUAL  TO  ZERO  -----  23 

1.  Analytical  Solution-  ---------------  23 

2.  Numerical  Solution  ----------------  24 

B.  INTERNAL  FORCING  EQUAL  TO  ZERO  ---------  -  27 

1.  Analytical  Solution-  ---------------  27 

2.  Numerical  Solution  ----------------  30 

C.  GENERAL  CASE _______  33 

D.  BAROCLINIC  CASES 35 

1.  Internal  Heating  Equal  to  Zero  __________  35 

a.  Positive  Temperature  Gradient-  --------  38 

b.  Negative  Temperature  Gradient-  --------  38 

2.  Influence  of  Internal  Heating-  ----------  43 

IV.  CONCLUSIONS 50 

LIST  OF  REFERENCES 52 

INITIAL  DISTRIBUTION  LIST 53 


LIST  OF  FIGURES 


1.  Two  dimensional  model  along  36°N  with  Az  =  40m  and 

Ax  =  80km 17 

2.  Profile  of  the  rate  of  conversion  of  available 
potential  energy  to  kinetic  energy „   From  Gill 

Green  and  Simmons  (1974)-  -----------------   20 

3.  Mean  velocity  profiles  from  Gill,  Green  and  Simmons 

(1974)  used  to  determine  temperature  profiles  -------  21 

4.  Profile  of  vertical  derivative  of  w'T*  showing  a 
large  upward  flux  in  the  first  500m.   Units  are 

°C  sec"1  x  10"7 22 

5.  Analytic  solution  to  the  case  in  which  wind  forcing 

is  equal  to  zero-  ---------------------   25 

6.  Initial  temperature  (°C)  profile  for  numerical  model-  -  -  -   26 

7.  Intermediate  numerical  solution  to  case  in  which 
wind  forcing  equal  to  zero,   a)  28  years  b)  55 

years   c)  82  years   d)  165  years-  -------------   28 

8.  Analytic  solution  to  the  case  in  which  heating  by  meso- 

scale  eddies  is  equal  to  zero  with  w  =  10  cm  day~l  -  -  -  -   31 

o 

9.  Note  the  0°  lines  (no  temperature  gradient)  which 
represent  the  barotropic  solution  to  the  balance 
between  w(8t/8z)  and  K(82T/8z2).   From  Stommel 

and  Webster  (1962) 32 

10.  Numerical  solution  to  the  case  in  which  heating 
by  mesoscale  eddies  is  equal  to  zero  with 

w  =  -10  cm  day-1  ---------------------  34 

o  J 

11.  Numerical  solution  to  the  general  case  with 

w  =   10  cm  day~l  ---------------------  36 

o  2 

12.  Numerical  solution  to  the  general  case  with 

w  =  -10  cm  day~l  ---------------------  37 

o  2 


13.  Initial  temperature  field  with  a  5°C  positive 

gradient  in  the  east-west  direction-  -----------   39 

14.  Numerical  solution  for  OT/9x)Q  >  0  and  w  = 

10  cm  day"1,  0/8z)  (w'T')  =0 ° 40 

15.  Numerically  calculated  internal  vertical  motion 

for  (8t/8x)  >  0  and  w  =  10  cm  day-1,   w  is  in 

1       o         o 
cm  sec"1  x  10"*-  ---------------------   41 

16.  Initial  temperature  field  with  a  5°C  negative 

gradient  in  the  east-west  direction-  -----------   42 

17.  Numerically  calculated  internal  vertical  motion 
in  cm  sec"1  x  10" ->  for  (dT/dx)    <   0  and  w  = 

10  cm  day"1-  -----------------------   44 

18.  Internal  vertical  motion  for  (8T/8x)  <  0  and 

—  1  —  1      —  r 

w  =  -10  cm  day   .   w  is  in  cm  sec   x  10" D-   -------   45 

o 

19.  Numerical  solution  for  QT/9x)  <  0  and  w  <  0. 

Temperature  is  in  °C.  0/8z)  (w'T')=  0----------   46 


20.   Internal  vertical  motion  for  (8T/8x)^<  0  and 

1 

o 


1  •      c'  O 

w     >   0.      w  is    in  cm   sec   ixl0-' 48 


21.      Numerical   solution  for    Ot/8x)q<    0  and  w     = 

10  cm  day"1-   -----------------------        49 


ACKNOWLEDGEMENTS 

The  author  wishes  to  express  her  thanks  to  Dr.  Robert  L.  Haney 
for  his  guidance  given  throughout  the  preparation  of  this  thesis. 

Special  appreciation  is  given  to  the  night-shift  personnel  of 
the  W.  R.  Church  Computer  Center  and  to  a  number  of  understanding 
friends  who  have  been  patient  with  me  during  the  writing  of  this 
thesis. 


I.   INTRODUCTION 

In  recent  years  scientists  have  found  that  highly  energetic,  rotat- 
ing water  masses  exist  in  the  ocean  below  the  mixed  layer.   Eddies,  as 
these  masses  have  been  called,  appear  to  be  similar  to  the  high  and  low 
pressure  systems  in  the  atmosphere.   Their  origin  and  the  extent  of  their 
influence  on  the  temperature  structure  of  the  ocean  is  not  fully  known. 
Meteorologists  as  well  as  oceanographers  are  showing  an  increasing  in- 
terest in  these  eddies  as  a  means  of  obtaining  a  better  understanding  of 
the  parameters  of  air-sea  interactions. 

In  1970  a  group  of  scientists  began  preparations  that  would  lead  to 
documentation  of  the  existence  of  eddies.   The  Mid-Ocean  Dynamics  Experi- 
ment (MODE)  was  a  multi-national  experiment  conducted  in  the  area  between 
Bermuda  and  Florida.   Through  extensive  measurements  of  currents,  tem- 
perature, salinity  and  pressure  the  scientists  discovered  the  existence 
of  two  types  of  eddies;  deep  water  eddies  and  highly  energetic  eddies  in 
the  upper  1500  meters  of  the  ocean. 

It  is  the  purpose  of  this  paper  to  study  the  possible  influence  of 
the  upper-layer  eddies  on  the  thermal  structure  of  the  ocean  following 
the  methods  used  in  previous  studies  of  the  oceanic  thermocline.   These 
thermocline  theories  are  all  based  on  the  assumptions  of  geostrophic  and 
hydrostatic  balance,  incompressibility,  neglect  of  salinity  and  some 
form  of  the  steady  state  first  law  of  thermodynamics  which  in  most 
general  form  is 

3T     3t     3t    32T 

Utt—  +v-r—  +w-r—  =  < (1) 

3x     3y     9z      2 


Table  I  gives  a  list  of  symbols  and  constants  used  in  this  study.   The 
general  balance  expressed  by  (1)  is  a  balance  between  horizontal  and 
vertical  advection  and  vertical  diffusion  of  heat.   Stommel  and  Webster 
(1962)  have  considered  the  simple  barotropic  case  which  balances  verti- 
cal advection  and  vertical  diffusion.   Robinson  and  Stommel  (1958) 
ignore  the  effects  of  east-west  temperature  advection  (term  one  of  (1) ) 
while  Welander  (1971)  considers  a  balance  of  all  the  terms  in  equation 
(1).   In  each  of  the  above  studies,  the  effect  of  mesoscale  heating  was 
ignored. 

The  plan  of  this  study  is  to  consider  the  simplest  model  in  which 


the  effect  of  heating  by  mesoscale  eddies  given  by  the  term,  -*■— (w'T1), 
can  be  added  to  (1) .   If  the  energy  source  for  the  eddies  is  available 


potential  energy,  then  w'T'  and  its  vertical  derivative  will  not  be 
zero.   Since  the  eddies  have  been  verified  by  the  MODE  data,  it  there- 
fore becomes  important  to  examine  the  effect  that  this  term  may  have  on 
solutions  to  the  thermocline  problem.   For  this  study,  the  correlation 


term,  ^— (w'T1)  was  specified  according  to  theoretical  results  of  Gill, 
dz 

Green  and  Simmons  (1974)  for  a  baroclinically  unstable  eddy  of  wave- 
length 322  km. 

Since  the  solution  to  the  heat  equation  used  for  this  study  depends 


on  the  amplitude  of  w'T',  its  profile  with  depth  and  on  the  wind  forced 
Ekman  vertical  velocity,  cases  were  run  to  test  each  of  these  factors. 
Analytical  and  numerical  solutions  are  found  for  the  cases  in  which  the 
Ekman  wind  and  then  the  mesoscale  heating  are  set  equal  to  zero  respec- 
tively. A  numerical  solution  is  found  when  both  the  wind  forcing  and 
the  internal  heating  are  included.   In  all  cases  horizontal  temperature 
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TABLE  I.   Symbols  and  Values  of  Constants 

t  time 

w  vertical  velocity 

T  perturbation  temperature 

R  radius  of  Earth 

<J>  latitude,  36°N 

X  longitude 

ft     angular  speed  of  Earth's  rotation,  2tt  day 

-4     -1 
a  thermal  expansion  coefficient,  2.0  x  10    (°K) 

f     coriolis  parameter,  2ft  sin  <J> 

v     northward  velocity  component 

u     eastward  velocity  component 

2    -1 
K     vertical  eddy  viscosity  and  conductivity,  .5  cm   sec 
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advection  is  ignored.   The  results  obtained  in  this  study  indicate 
that  the  effect  of  baroclinically  unstable  eddies  on  the  ocean  tempera- 
ture profile  is  significant. 
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II.   EQUATION  AND  MODEL  DEVELOPMENT 

A.   HEAT  EQUATION 

The  basic  heat  budget  equation  for  the  ocean  is 


§  ■  Q/ocp  (2> 


where  Q   represents  the  convergence  of  net  incoming  and  outgoing  radia- 
tion.  It  is  assumed  that  radiative  exchange  takes  place  in  the  first 
few  meters  of  the  ocean.   The  model  designed  for  this  study  is  such  that 
the  level  of  z  =  0  is  at  the  bottom  of  the  surface  Ekman  layer.   This  im- 
plies that  within  the  ocean  itself 


dT 

d^=°  (3> 


Expanding  (3)  and  writing  it  in  flux  form  gives 


oT 

||  +  V  •  C'.T)  =  0  (4) 


Let  W  =  '7  +  W   and  T  =  T  +  T'   where  at  this  point  the  primes  denote 
all  turbulent  motions  including  the  baroclinic  mesoscale  eddies.   A  time 
average  sufficient  to  remove  all  the  turbulent  motions  (say  120  days)  is 
performed  with  the  result 

8T 


—  +  \y  •  VT  +  V  •  Ov'T')  =  0  (5) 

dt 
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Expanding  this  into  component  parts  and  making  the  simplifying  assump- 
tion that  advection  and  heat  transfer  by  turbulence  in  the  north-south 
and  east-west  directions  can  be  ignored  gives  the  following  equation 

3T   -  8T 


3t  +  W-97  +  -3l  (W'T,)  =  °  (6) 


The  turbulent  vertical  transport  of  heat, (w'T'),  can  be  thought  of 
as  being  accomplished  by  two  distinct  physical  processes.   The  first  is 
through  transport  by  baroclinic  mesoscale  eddies  represented  by  the 


term  (w'T1)        ,   where  the  subscript  of  mesoscale  will  be  assumed 
mesoscale 

for  the  rest  of  this  paper.   The  second  means  of  heat  transfer  is  by  all 
other  turbulent  motions.   This  transfer  is  assumed  to  be  adequately  re- 
presented by  vertical  diffusion  of  heat  down  the  mean  temperature 

3T 
gradient;  -K  -~—  .   The  final  equation  to  be  solved  is  therefore 

0  Z 

3t   -  3t   3  ,—t-t-         9t\ 

lhT+w  87  +  al  (WT  "K  3l>  =  °  (7) 


where  w'T'  now  represents  the  heat  transport  by  mesoscale  eddies. 

B.   M  EQUATION 

This  study  follows  the  approach  of  Welander  (1971)  who  considered 
the  solution  to  a  similar  heat  equation  for  the  thermocline  region, 
but  without  the  heating  due  to  mesoscale  eddies.   The  steady  state, 
geostrophic,  hydrostatic  equations  in  a  region  away  from  side  bound- 


aries are 


-fv  =  - P,  ,  fu  =  -  —  P ,  ,  P  =  gaT 

Rcos(f)  A  R  <p  z 

—  cos*  {uA  +  (vcoscb)  ,}  +  w  =0 
R  (p     z 

u        v 

T,  +—  T^+wT   =KT 


RcoscJ)   A    R   (J)       z       zz 
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where  the  variable  P  and  T  are  perturbation  pressure  and  perturbation 
temperature.   Salinity  is  neglected  and  the  assumptions  of  a  spherical 
earth,  ocean  depth  H  small  compared  to  the  radius  of  the  earth  and  con- 
stant density  in  the  horizontal  direction  are  made.   Following  Needier 
(1967) ,  Welander  derived  a  non-linear  equation  for  pressure  of  fourth 
order  and  second  degree.   Noting  the  difficulty  in  solving  the  equation 
for  pressure  he  defined  a  new  function 


M(A,<}>,z)  =  /  pdz  +  2ftR  sin  <$>  J  w  dA 


(8) 


which  is  a  first  integral  of  pressure.   The  function  also  incorporates 
a  boundary  condition  since  w  is  the  vertical  velocity  at  the  bottom  of 
the  surface  Ekman  layer.   Welander  points  out  that  in  general  (u,v)  ^  0 
at  the  bottom  thus  implying  a  barotropic  component.   The  M  function  pro- 
vides the  following  definitions: 


u   ■  - — — — —  M  v   =  — — — — — ^—  M, 

2ftRsin<}>      <))z  2ftRsin<|>cos({)      Az 


w  =  = r—  M,    P=M        T  =  M 

„ri_2  .  2,   A        z  ga   zz 

2S2R  sin  c{> 


Substituting  these  definitions  into  the  heat  equation  (7)  developed 
for  consideration  in  this  thesis,  an  M  equation  similar  to  Welander  is 
obtained 

3  ...     .  l 


^•(M      )    + —  M.M  -KM  +   ga  -s-fw'T')    =    0  (9) 

dt     zz  ^„2    .    2,      A   zzz  zzzz  dz 

2S2R  sin  <p 


Letting 


m     =  e  do 

zz 
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and  assuming 

K  ?   K(z) 
there  results 


6.  +  = =-  M.6  -K0   +ga^(w'T')=0     (11) 

t    ^„2    2,   A  z     zz       dz 
2i3R  sin  <p 


The  first  term  of  (11)  is  proportional  to  the  local  time  change  of  tem- 
perature.  The  second  term  represents  the  vertical  advection  of  mean 
temperature  while  the  third  term  is  the  vertical  diffusion  of  the  mean 
temperature.   The  last  term  is  the  heating  due  to  mesoscale  eddies. 
Numerical  solutions  to  equations  (10)  and  (11)  for  a  steady  state  were 
obtained  using  the  initial  value  approach.   Analytic  solutions  were  also 
obtained  in  the  simplest  cases.   Before  discussing  analytic  and  numerical 
solutions  to  (10)  and  (11) ,  a  detailed  description  of  the  boundary  condi- 
tions will  be  given. 

C.   BOUNDARY  CONDITIONS 

The  physical  model  chosen  to  test  the  influence  of  ocean  eddies  is 
shown  in  Figure  1.   The  model  is  tested  in  a  region  of  the  Pacific  Ocean 
away  from  any  strong  boundary  currents,  along  36°N.   The  domain  was 
chosen  to  be  large  compared  to  the  size  of  an  eddy  and  to  be  sufficiently 
small  so  that  the  atmospheric  forcing  could  be  considered  approximately 
uniform.   A  flat  bottom  was  assumed  for  the  two  dimensional  model.   For 
initial  test  cases  the  boundary  conditions  were  independent  of  longitude 
making  this  a  one-dimensional  model  (solution  independent  of  x) . 

Since  equation  (9)  is  of  fourth  order,  four  boundary  conditions  can 
be  applied  to  the  problem.   An  advantage  of  using  Welander's  M  function 
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*    X  ,U    (east) 


•Figure   1.      Two  dimensional  model   along   36°N  with  Az  =   40m  and 
Ax  =   80m. 
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is  that  part  of  the  boundary  conditions  are  incorporated  in  the  function 
itself.   In  practice  it  was  convenient  to  solve  (10)  and  (11)  as  the 
system 


MAzz  "  9A  (12) 


ae 


M,9   -Kfl,  +  get  ■£-  (w'T1)  =  0  (13) 


3t   on„2   2.      A  z     zz   ^   8z 
2S2R  sin  <J> 


These  are  solved  conveniently  with  boundary  conditions  on  M,  at  the 
ocean  top  and  bottom  in  equation  (12)  and  on  0  at  the  ocean  top  and 
bottom  in  equation  (13) . 

The  model  has  a  rigid  lid  with  w  prescribed  at  the  top  and  bottom. 

2    2 
At  z  =  0,  M,   at  the  surface  becomes   2ftR  sin  (b  w   where  w   is  a  pre- 
A  o  o 

scribed  vertical  motion  produced  by  the  wind  at  the  bottom  of  the 

surface  Ekman  layer.   At  z  =  -H  ,  the  boundary  condition  is  w  =  0  , 

which  requires   M,  =  0  .   The  thermal  conditions  at  the  top  and  at  the 

bottom  are  to  prescribe  the  temperatures   T   and  T    respectively. 

Thus,  at  the  surface   0  =  get  T   and  at  the  bottom,  0  =  ga  T  „  . 

o  -H 

D.   EDDY  HEAT  TRANSFER 

In  this  thesis  no  attempt  has  been  made  to  parameterize  the  eddy 


transfer  of  heat,  w'T1,  in  terms  of  the  dependent  variables.   No  direct 
measurements  of  the  heat  transfer  were  available  so  it  was  necessary  to 
resort  to  theory  to  obtain  an  estimate.   In  a  paper  by  Gill,  Green  and 
Simmons  (1974)  a  theoretically  determined  rate  of  conversion  of  avail- 
able potential  energy  to  eddy  kinetic  energy  is  developed  as  a  function 


of  z  in  a  vertically  continuous  model.   The  values  of  w'T'  presented  by 
Gill  et  al  are  those  that  would  exist  if  the  mesoscale  eddies  were  the 
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result  of  in  situ  baroclinic  instability.   As  such,  the  w'T*  which  is 
used  in  this  paper  is  based  on  a  given  horizontal  and  vertical  profile 


of  temperature.   Figure  2  from  Gill  e_t  a_l  shows  w*T*  as  a  function  of 
depth  for  eddies  of  four  different  wavelengths.   The  amplitudes  asso- 
ciated with  these  profiles  are  based  on  the  assumption  that  all  of  the 
available  potential  energy  in  the  mean  flow  is  converted  into  the 
eddies.   This  study  can  be  thought  of  as  a  test  of  this  assumption  as 
well  as  a  test  of  the  in  situ  baroclinic  instability  theory.   Figure  3 
shows  the  mean  velocity  which  determines  the  temperature  profile. 

For  analytic  treatment  in  this  thesis  the  heat  flux  due  to  meso- 
scale  eddies  was  expressed  as  an  exponential  function. 


z/D 


w'T'  =  A  e  (14) 

where  A  ,  .0357  °C  cm  sec   ,  is  the  amplitude  and  D  ,  100  m  ,  is  the 


characteristic  depth  of  energy  transfer.   The  profile  of  ■*—  (w'T1)  is 
shown  in  Figure  4.   This  corresponds  in  amplitude  and  profile  to  that 
shown  in  (b)  of  Figure  2.   The  large  amplitude  indicates  an  extremely 
strong  upward  flux  of  heat  in  the  upper  500  m  of  the  ocean. 
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Fig.  12.  The  rate  of  conversion  ol available  potential  enenty  to  eddy  energy,  shown  n*  a  function 
of  depth.  The  scale  is  chosen  so  that  the  integral  with  respect  to  depth  on  m)  gives  ii  in  W  m_;  $-. 
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Figure  2.  '  Profile  of  the  rate  of  conversion  of  available  potential 
energy  to  kinetic  energy.   From  Gill,  Green  and  Simmons 
(1974). 
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Fig.  ".  Mean  velocity  profiles  C(z)  for  which  a  stability  analysis  has  been  carried  out.  For  clarity, 
profiles  [4j  and  [51  have  been  displaced  to  the  right  and  profile  [  1 J  is  drawn  twice.  Note  that  the  scale 
of.-  is  iimitec  :o  leys  :nan  half  the  total  depth.  The  mean  currents  below  2000  m  are  very  small. 


Figure  3.   Mean  velocity  profiles  from  Gill,  Green  and  Simmons 
(1974)  used  to  determine  temperature  profiles. 
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III.   ANALYTIC  AND  NUMERICAL  RESULTS 

The  solutions  to  equations  (12)  and  (13)  are  dependent  on  the  inter- 
nal heating  and  the  surface  boundary  conditions  of  temperature  and  verti- 
cal velocity.   In  order  to  aid  in  the  interpretation  of  the  more  general 
numerical  solutions,  analytical  solutions  were  obtained  for  a  variety  of 
the  simpler  cases  first. 

The  first  case  to  be  examined  was  that  of  no  Ekman  forcing  and  the 
next  case  was  that  of  no  internal  mesoscale  eddy  heating.   The  final  case 
was  the  most  general  case  which  included  both  Ekman  forcing  and  heating 
by  mesoscale  eddies.   Both  upward  and  downward  vertical  velocities  were 
tested  where  applicable. 

A.   VERTICAL  MOTION  AT  THE  SURFACE  EQUAL  TO  ZERO 
1.   Analytical  Solution 

In  the  first  case  the  vertical  motion  at  the  bottom  of  the  surface 
Ekman  layer,  w  ,  was  set  equal  to  zero.   If  the  surface  temperature,  T  , 
is  independent  of  x,  the  solution,  0,  is  independent  of  x.   Therefore  by 
equation  (12) ,  M,  =  0  everywhere.   The  steady  state  form  of  equation  (13) 
thus  becomes  (using  (14)) 

92T    3  — r—r,         A   z/D  ,,_, 

K  =  -*— (w'T')  =  —  e  (15) 


dz 


3z         D 


This  is  the  equation  for  a  local,  one-dimensional  model  without  vertical 
advection.   The  analytic  solution  to  (15)  was  obtained  by  integrating  the 
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diffusion  term  twice  and  applying  the  boundary  conditions  of  T(0)  =  Tq 
and  T(-H)  =  T    . 


T(z)  -  T  +  |  {(T  -  T   »  -  ^(l-e-H/D)}  -  ^(l-e2/D)   (16) 
O    rl     O    ~ri      IC  K 


The  solution  (16)  is  shown  in  Figure  5.   The  temperature  is 
forced  to  a  value  of  -554°C  at  approximately  360  m  below  the  surface. 
The  extremely  unrealistic  nature  of  this  solution  indicates  that  if  the 
open  ocean  contains  mesoscale  eddies  with  the  above  heating  profile, 

then  the  thermal  balance  cannot  be  maintained  by  simple  downward  heat 

2    -1 
diffusion  with  K  =  . 5  cm  sec 

2.   Numerical  Solution 

In  order  to  test  the  numerical  model  which  will  be  used  later 

for  the  more  general  case,  a  numerical  solution  to 


T  =  KT   -  t-Iw'T')  (17) 

t     zz    dz 


was  obtained.   The  numerical  method  is  the  initial  value  integration 
starting  from  an  initial  temperature  field  given  by 

T(z)  =  TieZ/dl  -  T2eZ/d2  (18) 

The  initial  profile  was  defined  to  meet  the  stable  condition  of 

8T 

-r—  _>  0.   For  this  to  hold  the  following  constants  were  chosen;  T  = 

22.5°C,  T  =  2.5°C,  d  =  900  m  and  d  =  100  m.   Figure  6  shows  this 

profile. 

The  general  numerical  finite  difference  scheme  to  solve  (17) 

used  a  centered  differencing  in  space  and  time  except  for  the  diffusion 
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Figure  5.   Analytic  solution  to  the  case  in  which  wind  forcing  is 
equal  to  zero. 
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Figure  6.   Initial  temperature  (°C)  profile  for  numerical  model, 
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term  which  utilized  a  forward  time  step.   The  finite  difference  analogue 


is 


Tn+1  =  Tn_1  +  2At(KTn  1  -   ga  JLcw'T'))  (19) 

ZZ  02. 


where  At  =  10  days,  Az  =  40  m,  Ax  =  80  km.   The  final  solution  was 
reached  after  integrating  for  approximately  2000  years  and  is  in  agree- 
ment with  the  analytic  solution.   Intermediate  solutions  showing  the 
approach  to  steady  state  are  shown  in  Figure  7.   It  can  be  noted  that 


the  forcing  term,  -*—  (w'T')f  immediately  works  to  cool  the  temperature 
in  the  top  100  meters  and  then  the  diffusion  term  acts  to  redistribute 
the  cooler  water  throughout  the  entire  layer.   Solutions  to  the  baro- 
clinic  case  forced  by  a  thermal  gradient  are  considered  in  section  D 
below. 

B.   INTERNAL  FORCING  EQUAL  TO  ZERO 
1.   Analytical  Solution 

In  the  second  case  to  be  tested  the  heating  due  to  mesoscale 


eddies,  w'T',  was  set  equal  to  zero.   This  required  that  there  be  a 
balance  between  heat  diffusion  and  vertical  temperature  advection  for 
steady  state.   From  (12)  and  (13)  the  equations  to  be  solved  were 

2    2 

2ftR  sin  (fr  m                             ___ 

*■  w   -  T,                            (20) 

ga      zz  X 

-  32T     3T  ,__. 

-  W  -5-  (21) 


9z 


8z 


If  the  surface  conditions,  T  and  w  ,  are  independent  of  x,  equation  (2) 

o      o 

reduces  to  w   =0  with  solution 
zz 

w  =  wq(1  +  |)  (22) 
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Figure  7.   Intermediate  numerical  solution  to  case  in  which 
wind  forcing  equal  to  zero,   a)  28  years  b)  55 
years   c)  82  years  d)  165  years. 
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The  equation  to  solve  is  thus 


82T 


Z.  9T 


9z 


=  w  (i  +  ^)  — - 
O      H   dz 


(23) 


with  boundary  conditions  T(0)  =  T  and  T(-H)  =  T  „ 

o  -H 


This  is  the  local, 
-1 


one-dimensional  equation  with  Ekman  forced  vertical  motion  of  10  cm  day 


but  no  heating  due  to  the  mesoscale  eddies. 

9T 
Let  y  =  -tt—   ,  then  (23)  becomes 
dz 


y 


w 

-*  <1  + 1) 

K        H 


Thus, 


w        2 
£ny  =  f  (Z  +  3?  +  °1  ' 


or 


3t         ir(z  +  z2/2H) 

y  =  U  =  C2  e 


Integrating  from  -H  to  Z; 


w 


T(Z)  -  T_H  =  c2 


/.*•  * 


n2/2H) 


dn 


-H 


where  n,   is  the  variable  of  integration.   Evaluating  c  by  requiring 
T(0)  =  T  ,  gives  the  final  solution 


T(Z)  =  T  „  + 


T  _  T 
O     -H 


O   w 


u e 


— (n+n  /2H) 


dn 


w 

5  -?(n+n  /2H) 


K 


dn 


-H 


(24) 


The  solution  is  not  in  closed  form  but  is  given  by  an  integral 
which  must  be  calculated  numerically.   The  results  of  equation  (24) 
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being  obtained  by  numerical  integration  using  the  trapezoidal  rule  are 
shown  in  Figure  8.   Notice  that  the  stratification  is  confined  to  the 
upper  500  m. 

Stomrael  and  Webster  (1962)  have  analyzed  a  similar  case  in 
which  vertical  heat  diffusion  and  heat  advection  are  balanced.   Their 
solution  in  the  barotropic  case  shows  upward  vertical  velocity  a  linear 
function  of  depth,  as  in  (22) ,  but  the  temperature  is  isothermal.   The 
similarity  transform  used  to  solve  their  coupled  equations  does  not 
allow  vertical  temperature  stratification  in  the  barotropic  case. 
Figure  9  shows  the  case  from  Stommel  and  Webster  (1962)  in  which  the  w 
and  T  lines  denoted  by  0°  (implying  a  zero  horizontal  temperature  gradi- 
ent) are  for  the  barotropic  case.   The  stratified  temperature  profile 
obtained  in  this  study  is  allowed  because  equations  (12)  and  (13)  are 
solved  without  any  restrictions  imposed  by  a  similarity  transformation. 
2.   Numerical  Solutions 

In  order  to  check  the  numerical  model,  the  same  case  of  no  heat- 
ing by  mesoscale  eddies  was  solved  numerically.   By  assuming  that  the 

8T 
local  temperature  change  with  time,  -tt—  ,  is  not  equal  to  zero,  an  initial 

d  t 

value  time  integration  was  run  on  the  computer.   The  conditions  are 
identical  with  those  leading  to  equation  (23)  so  the  equation  to  inte- 
grate was 


8T     92T      ,.    Z,  3T  ,„_, 

3zJ 


-irr  -    r*   -  w  \±    -r   —j    -~— 

dt      ^2  O      H   dz 


The  finite  difference  analogue  to  this  was 


mn+l   n-1   _.  .   n-1      ,,   Z.    nN    ,_„. 

T    =  T    +  2At(KT     -  w  (1+  — )  T   )    (26) 

zz     o    H    z 
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Figure  8.  Analytic  solution  to  the  case  in  which  heating  by  meso- 

scale  eddies  is  equal  to  zero  with  w  =  10  cm  day-1. 
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The  model  was  run  with  a  time  step  of  ten  days  for  1000  years.   A  steady 
state  was  reached  after  500  years  for  the  case  when  w  >  0  (upward 
motion) .   The  numerical  solution  agrees  with  the  analytically  obtained 
solution  shown  in  Figure  8.   Since  the  numerically  and  analytically  ob- 
tained solutions  agree  in  these  simple  cases,  it  is  reasonable  to  trust 
the  accuracy  of  the  numerical  solution  in  the  more  general  cases  con- 
sidered below. 

A  numerical  solution  to  the  case  with  no  internal  heating  was 
run  with  downward  vertical  velocity  of  the  same  magnitude  as  used  before; 
-10  cm  day   .   Figure  10  shows  an  isothermal  condition  to  great  depth 
with  a  strong,  deep  thermocline.   It  appears  that  the  vertical  tempera- 
ture advection  term  works  in  the  same  sense  as  the  vertical  diffusion. 
Such  a  profile  is  clearly  not  observed  in  the  real  ocean. 

C.   GENERAL  CASE 

When  both  physical  effects  of  heating  by  mesoscale  eddies  and  wind 

forcing  are  included  recall  that  the  equations  to  solve  are  (12)  and 

(13) .   In  the  case  in  which  T  and  w  are  independent  of  x,  the  solution 

o      o 

is  independent  of  x,  and  so  w  is  given  by  (22) .   Equation  (13)  becomes, 
with  substitution  of  (22) , 

86         2.  A  3 


xr  +  wll  +  3  fl  -  k6   +  get  -£-  (w»T')  =  0  (27) 

dt  o     Hz      zz       dz 

The  steady  state  version  of  this  equation  is  similar  to  (23)  with  the 
addition  of  the  non-homogeneous  term  and  an  analytic  solution  is 
possible.   Due  to  the  complexity  of  the  anticipated  solution  to  equation 
(27)  an  analytic  solution  incorporating  all  three  terms  was  not  obtained. 
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Figure  10.  Numerical  solution  to  the  case  in  which  heating  by_meso- 
scale  eddies  is  equal  to  zero  with  w  =  -10  cm  day   . 
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Only  the  initial  value  numerical  method  of  solution  was  considered. 

Using  a  time  step  of  ten  days  steady  state  solutions,  shown  in  Figures 

11  and  12,  were  reached  in  approximately  500  years  for  the  cases  of 

w  >  0  and  w  <  0  respectively, 
o  o 

When  w  is  upward  (Figure  11)  the  initial  temperature  change  due  to 

the  forcing  function  is  similar  to  the  case  when  w  =  0  (Figures  5  and 

7) .   The  first  250  m  are  cooled  rapidly  but  then  the  influence  of  the 

vertical  advection  slows  the  cooling  which  is  then  confined  to  the  upper 

region.   The  profile  is  more  realistic  than  that  given  in  Figure  5  for 

the  case  w  =  0;  however  the  temperature  is  still  an  unrealistic  -174 °C 

at  200  m. 

When  w   is  downward,  the  result  is  seen  in  Figure  12.   Again  cooling 
o 

takes  place  quite  rapidly  due  to  the  influence  of  the  mesoscale  eddies. 
The  downward  vertical  velocity  works  with  the  vertical  diffusion  so  the 
deeper  waters  also  experience  cooling.   It  has  been  shown  that  none  of 
the  local  solutions  which  contain  heating  by  mesoscale  eddies  are  realis- 
tic.  Further  conclusions  will  be  reached  after  examining  the  baroclinic 
cases  below. 

D.   BAROCLINIC  CASES 

1.   Internal  Heating  Equal  to  Zero 

The  final  series  of  test  cases  were  designed  to  determine  the  in- 
fluence of  heating  by  mesoscale  eddies  in  a  baroclinic  ocean.  This  test- 
ing was  accomplished  by  comparing  the  results  obtained  with  and  without 


the  mesoscale  eddy  heat  transfer  term,  ■*—  (w'T1).   The  initial  value 

oz 

approach  to  the  solution  was  used  with  a  temperature  field  having  a  five 
degree  surface  gradient. 
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Figure  11.   Numerical  solution  to  the  general  case  with 

w  =  +10  cm  day"1. 
o  J 
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Figure  12.  Numerical  solution  to  the  general  case  with 

w  =  -10  cm  day-1. 
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The  equation. 


T(x,z)  =  (T  +  5.(x/L))eZ/dl  +  ^^^2  (28) 


where  L  is  the  east-west  extent  of  the  model  (2000  km) ,  was  used  to 
generate  an  initial  field  with  either  a  positive  or  a  negative  east- 
west  gradient.   The  values  of  T(x,0)  and  T(x,-H)  obtained  from  (28) 
constituted  the  fixed  thermal  boundary  conditions.   With  the  addition 
of  an  east-west  temperature  gradient,  the  internal  vertical  motion  is 
now  a  function  of  depth  and  temperature.   Because  equation  (13)  is  now 
non-linear,  M,  must  be  calculated  at  each  time  step. 

a.  Positive  Temperature  Gradient 

Figure  13  shows  the  initial  temperature  field  with  a  posi- 
tive gradient;  warmer  water  to  the  east.   The  numerical  solution  to 
equations  (12)  and  (13)  is  shown  in  Figure  14  for  the  case  when  the  sur- 
face Ekman  forcing  is  greater  than  zero,  w  =  10  cm  day   .   An  unrealis- 
tically  deep  isothermal  layer  is  depicted.   This  is  also  the  case  for 
w   less  than  zero.   A  strong  interior  downward  vertical  motion,  shown  in 
Figure  15,  is  created  by  the  temperature  gradient  in  both  of  the  above 
cases.   This  sinking  motion  works  with  the  diffusion  to  cause  a  packing 
of  the  isotherms  near  the  bottom  of  the  ocean.   With  no  other  forces  act- 
ing to  counter  the  double  influence  of  the  downward  heat  advection  and 
diffusion,  the  resolution  of  the  model  is  insufficient  to  handle  the 
strong  gradients  that  finally  develop. 

b.  Negative  Temperature  Gradient 

The  initial  field  with  a  five  degree  negative  east-west 
gradient  at  the  surface  is  shown  in  Figure  16.   When  the  model  was  run 
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Figure  13.   Initial  temperature  field  with  a  5°C  positive  gradient 
in  the  east-west  direction.  . 
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Figure  14.   Numerical  solution  for  (3t/3x)  >  0  and  w  =  10  cm  day"1, 
(3/3z) (w'T')  =  0.   Temperature  is  in  °C.° 
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Figure  16.   Initial  temperature  field  with  a  5°C  negative  gradient 
in  the  east-west  direction. • 
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with  cyclonic  wind  conditions  generating  an  Ekman  layer  upward  motion 
of  w  =  10  cm  day   ,  the  resulting  interior  vertical  motion  is  strongly 
positive  throughout  the  model  (Figure  17) „   The  resulting  vertical  ad- 
vection  is  stronger  than  the  downward  diffusion  and  causes  a  packing  of 
the  isotherms  near  the  surface. 

When  anticyclonic  wind  conditions,  generating  downward  ver- 
tical motion  at  the  base  of  the  surface  Ekman  layer,  are  applied  the 
interior  vertical  motion  becomes  positive  only  below  250  m.   Its  maxi- 
mum value  is  smaller  than  in  the  previous  case  in  which  w  was  positive. 
Figure  18  shows  a  profile  of  the  vertical  motion  taken  from  the  central 
region  of  the  model.   For  this  case  the  vertical  temperature  advection 
does  not  dominate  the  downward  diffusion  and  only  a  slight  packing  of 
isotherms  occurs.   A  steady  state  solution,  reached  after  50  years,  is 
shown  in  Figure  19  with  a  well  defined  thermocline  present.   This  is 
clearly  representative  of  the  typical  situation  below  an  atmospheric 
subtropical  high  pressure  system  with  near  surface  temperatures  gradual- 
ly decreasing  eastward  from  the  center  of  the  oceanic  anticyclonic  gyre. 

Without  internal  heating  the  effect  of  a  temperature 
gradient  is  as  follows.   When  the  east-west  surface  temperature  gradient 
is  positive,  the  internal  vertical  motion  becomes  negative  regardless  of 
the  surface  value  for  w  .   The  resulting  temperature  profiles  are  similar 
to  Figure  10.   A  negative  surface  gradient  produces  positive  interior 
values  for  w  and  the  baroclinic  solutions  without  internal  heating  are 
similar  to  Figure  8. 

2.   Influence  of  Internal  Heating 

The  influence  of  heating  by  mesoscale  eddies  was  tested  for  the 
baroclinic  case  with  a  negative  horizontal  temperature  gradient  and 
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Figure  18.   Internal  vertical  motion  for  (3t/3x)0<  0  and  w 
-10  cm  day"  .   w  is  in  cm  sec~l  x  10"^. 
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Figure  19.      Numerical   solution  for    (3T/3x)ri<   0   and  wq<    0, 
Temperature   is   in   °C.       (3/3z)     (w'T')=   0. 
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positive  Ekman  forcing.   Equations  (12)  and  (13)  were  solved  using  the 

initial  value  approach  with  the  interior  vertical  motion  calculated  at 

each  time  step  as  a  function  of  depth  and  temperature.   After  40  years 

of  integration,  using  a  At  of  ten  days,  the  result  was  a  vertical  motion 

profile  shown  in  Figure  20.   Recall  that  in  the  local  solution  case 

which  included  the  eddy  term  with  w  >  0  the  interior  value  of  w  was 

o 

simply  a  linear  function  of  depth  (Section  III,  C) . 

The  accompanying  solution  for  temperature  in  the  baroclinic  case 

is  shown  in  Figure  21.   In  this  case  the  influence  of  the  eddies  was  to 

act  in  the  same  sense  as  the  interior  vertical  motion  and  cool  the  near 

surface  layer  to  an  unrealistic  -142°C.   This  is  comparable  to  the  local 

case  with  w  >  0  (Figure  11) ,  in  which  the  minimum  temperature  reached 
o 

was  -174 °C.   It  is  evident  that  the  baroclinicity  acts  to  moderate  the 
cooling  influence  of  the  mesoscale  eddies  but  only  to  a  small  extent. 


Gill  et  al^  showed  that  the  amplitude  used  for  w'T*  in  this  paper  was  re- 
quired in  order  that  in  situ  baroclinic  instability  produce  eddies  of 
the  intensity  observed  in  the  open  ocean.   It  is  of  interest  to  determine 


the  sensitivity  of  the  solution  to  the  amplitude  of  w'T'.   Test  cases 


were  run  in  which  the  amplitude  of  w'T'  was  varied.   It  was  found  that 

reducing  the  amplitude  by  two  orders  of  magnitude  was  required  to  produce 

a  mean  temperature  profile  close  to  that  observed  in  the  real  ocean. 

Cases  were  also  run  with  internal  heating  and  the  following  con- 

9T  8T  8T 

ditions;  t~-  <  °  with  w  <  0,  -r— —  >  0  with  w  >  0  and  -r— —  >  0  with  w  < 
dx  o      dx  o         dx  o 

0.   In  all  of  these  cases  the  resolution  of  the  model  was  insufficient 
to  handle  the  strong  vertical  temperature  gradients  which  developed  near 
the  sea  surface.   Finer  vertical  resolution  and  thus  smaller  time  steps 
would  be  required  to  obtain  a  numerical  solution  for  any  of  the  above 
three  cases. 
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Figure   21.      Numerical   solution   for    OT/3x)Q  <   0  and  wq  - 


10  cm  day-1. 


49 


IV.   CONCLUSIONS 


It  is  obvious  that  the  real  ocean  does  not  freeze  as  a  result  of 
vertical  heat  transport  by  mesoscale  eddies.   There  are  a  number  of 
possible  explanations  for  the  excessive  cooling  observed  in  some  of  the 
test  cases  in  this  study.   In  deriving  a  model  that  was  simple  enough 
to  test  the  influence  of  eddies,  it  was  necessary  to  make  numerous 
assumptions.   With  the  motion  independent  of  y  and  with  the  temperature 

specified  at  the  top  and  bottom  of  the  model,  the  mesoscale  eddy  heat 

g  

transport  could  only  be  distributed  vertically  by  the  term  -r—  (w'T'). 

Were  there  an  advective  motion  in  the  meridional  direction,  the  north- 


ward eddy  heat  transport,  v*T',  could  have  been  distributed  by  the 


term  -s—  (v'T')  such  as  to  produce  much  less  drastic  effects  and  more 

3y 

realistic  temperature  profiles.   It  would  be  of  interest  to  include 
horizontal  temperature  advection,  v  -r—  ,  and  -r—  (v'T1)  and  then  run  the 
model  testing  the  sensitivity  of  those  terms. 

Since  the  mesoscale  heating  by  transient  eddies  was  determined 
from  an  analytic  consideration  of  energy  transfer  rate  by  Gill  e_t  al 
(1974) ,  it  is  possible  that  the  amplitude  of  that  transfer  was  overesti- 
mated.  As  stated  before,  Gill  e_t  al  are  estimating  a  maximum  eddy 
strength  by  assuming  that  all  of  the  mean  flow  available  potential 
energy  is  converted  into  eddy  energy.   Since  the  earth's  rotation  acts 
as  a  constraint  to  energy  conversion  and  since  the  most  favorable  condi- 
tions for  eddy  growth  are  found  where  the  isopycnals  slope  upward 
towards  the  equator,  i.e.,  colder  water  near  the  equator,  a  complete 
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conversion  of  mean  available  potential  energy  to  eddy  kinetic  energy 


may  not  be  possible.   Using  the  amplitude  for  w'T'  from  Gill  et.  al/ 
the  unrealistic  results  from  the  model  strongly  suggest  that  the  eddies 
observed  in  the  open  ocean  are  not  the  result  of  iji  situ  baroclinic 
instability. 

The  results  of  this  thesis  give  a  preliminary  indication  of  the 
potential  magnitude  of  influence  that  energetic  mesoscale  eddies  have 
on  the  temperature  structure  of  the  ocean.   It  is  evident  that  the 
heat  equation  in  the  ocean  should  include  mesoscale  eddy  heating  terms 
if  accurate  and  complete  solutions  are  to  be  obtained. 
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